Experimental evidence of solitary wave interaction in Hertzian chains 
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We study experimentally the interaction between two solitary waves that approach one to another 
in a linear chain of spheres interacting via the Hertz potential. When these counter propagating 
waves collide, they cross each other and a phase shift respect to the noninteracting waves is intro- 
duced, as a result of the nonlinear interaction potential. This observation is well reproduced by 
our numerical simulations and it is shown to be independent of viscoelastic dissipation at the beads 
contact. In addition, when the collision of equal amplitude and synchronized counter propagating 
waves takes place, we observe that two secondary solitary waves emerge from the interacting re- 
gion. The amplitude of secondary solitary waves is proportional to the amplitude of incident waves. 
However, secondary solitary waves are stronger when the collision occurs at the middle contact in 
chains with even number of beads. Although numerical simulations correctly predict the existence 
of these waves, experiments show that their respective amplitude are significantly larger than pre- 
dicted. We attribute this discrepancy to the rolling friction at the beads contacts during solitary 
wave propagation. 

PACS numbers: 43.25.-|-y 45.70.-n 
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INTRODUCTION 



As it was shown in the pioneering work of V. 
Nesterenko [H-Q, wave propagation in one dimensional 
assembly of beads has a rich content in physics due to 
the nonlinear nature of the interaction between spheres. 
If the alignment does not support any static load, lin- 
ear acoustic modes of propagation (phonons) are forbid- 
den and any infinitesimal small perturbation propagates 
as solitary waves. This fully nonlinear regime corre- 
sponds to a sonic vacuum Several studies have re- 
ported both experimentally and numerically a large num- 
ber of phenomena associated with this type of waves in 
Hertzian linear chains, ranging from solitary wave inter- 
action with boundaries [^-ul and wave mitigation in ta- 
pered chains IsMllll. t o energy localization in chains with 
mass defects p^ - [l4| and modification of wave propaga- 
tion in lubricated granular chains 15, 1^. The weakly 
nonlinear regime, for which the amplitude of a stress per- 
turbation is smaller but comparable to the imposed static 
stress, also exhibits well known characteristics. For in- 
stance, it has been shown that the dynamics of chain of 
weakly loaded spheres obey the Korteweg- de Vries equa- 
tion [1, 17, l3l, which is known in turn to support soli- 
tons, kinks and breathers as exact solutions Il9|. One of 



the remarkable features of the weakly nonlinear KdV soli- 
tons resides in their ability to interact nonlinearly with 
each other and escape a collision by conserving not only 
their momentum and energy, but also the full charac- 
teristics of the initial waves (shape, width, amplitude). 
The only trace of the nonlinear nature of the interac- 
tion is a slight post-collision phase shift of each wave, 
compared to the non- interacting configuration 17, [l8| . 



One of the main features of solitary wave is the nonlin- 
ear character of their interaction for which the linear su- 
perposition principle does not hold. Here, we investigate 
the interaction of two identical solitary waves traveling 
in opposite directions in a linear chain of spherical elas- 
tic beads. In the nonlinear work frame, the amplitude 
of the total wave during the interaction is larger than 
the simple sum of two individual wave. Since the wave 
speed is a nonlinear function of wave amplitude, when 
the interaction ends, both waves have covered the dis- 
tance where they interacted in a shorter time [lij . The 
later implies the existence of a phase-shift in the form of 
a time delay, as for KdV solitons in the weakly nonlin- 
ear regime. The observation of such a delay thus may 
serve to probe the exact nature of the strongly nonlinear 
solitary waves that propagate in one-dimensional weakly 
loaded granular structures. It has also been discussed 
that energy partition in chains - for counter propagating 
waves initiated at the borders at equal time and equal 
amplitude - with a even number of beads (where waves 
collide at the beads contact), differs from the case of a 
chain with odd number of beads (where the bead at the 
center does not move during the wave interaction). In the 
first case, the interaction generates an excess of potential 
energy at the bead's contacts which later contributes to 
the formation of secondary solitary waves (SSW) 0, 21 



The study of frontal interaction between solitary waves 
of significantly different amplitudes and of over-taking 
solitary waves traveling in the same direction both 
require setups having an important number of beads. In 
the present work, only frontal collision between waves of 
similar amplitude were considered, which has the advan- 
tage of dealing with short chains and accordingly small 
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dissipation. This paper is organized as follows: a brief 
description of the soliton-like solution for a one dimen- 
sional granular chain is first reminded, a full description 
of our experimental setup and observations is then given, 
and an analysis of experimental and numerical results 
are finally discussed in terms of solitary wave interaction 
time. 



ANALYSIS AND NUMERICAL SIMULATIONS 

The physical behavior of a chain of N beads un- 
der elastic deformation is described by using the 
pairs Hertz potential 23|, Uh = (2/5)k(5^/^ where 
S is the deformation of two consecutive spheres, 
K-i = {0 + 0'){R-^+R'-^y/^,9 = 3{l-iy^)/4Y and i? 



and R' are the radii of curvature at the contact. Y and v 
are the Young modulus and Poisson's ratio, respectively. 
Since the force felt at the contact is the derivative of 
the potential energy with respect to 6 {Fh = k6^^^ , 
and Fh = if the beads lose contact), it is possible 
to describe the dynamics of the chain by the following 
nonlinear coupled differential equations: 
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Here, time derivatives are denoted by overdots and m is 
the mass of each bead. Under the long-wavelength ap- 
proximation A ^ i? (where A is the characteristic length 
of the perturbation) for the case of vanishing external 
static load, one can obtain the continuum limit of Eq. ([T]) 
by replacing the discrete function Un±i{t) by the cor- 
responding Taylor expansion of the continuous function 
u{x±R,t). Keeping terms up to the fourth order in spa- 
tial derivatives, it is possible to obtain an equation for 
the strain ip — —dxU 



V'^/2 + (2/5)i?V'/'a..(V^/^)l (2) 



which admits an exact periodic solution in the form of 
a traveling wave, — x ~ vt), with speed v. This 
solution is: 
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The second term expression for the wave veloc- 
ity is obtained by relating the maximum force F^ 
to the maximum bead velocity V„i. Experimen- 
tally, the maximum force Fm, the time of flight T 
and the duration r are obtained here from fitting 
the loading part {—T<t — T<T) of the measured 
force as a function of time to Nesterenko's solution 
F{t) ^ Fm cos6 [(t - T)/t] (with 2t = 2RVW/v). 

As a way to complete the experimental information 
shown in the following sections, we perform numerical 
simulations by solving Eq. ([T]) for a chain of iV = 25 
(or N = 26) beads, using a fourth order Runge-Kutta 
algorithm, with a time step smaller than the minimum 
physical time in the system and which fulfills the energy 
conservation up to relative variations of order 10~^. In 
our numerical simulations we have taken viscoelastic 
dissipation into account by using a nonlinear viscoelastic 
solid model ^ (F„ = {3/2)r]K6'^/^6, with viscous relax- 
ation time 7] — 0.58/is, derived for this set of beads, from 
the measurement of restitution coefficient Q). To initi- 
ate the waves, initial velocities are applied over the first 
and last beads {+Vo and —Vq, correspondingly). The re- 
sult of our simulations allow us to estimate the force felt 
by all contacts as function of time and provide a direct 
comparison with the forces measured in our experiments. 



EXPERIMENTAL SETUP 
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FIG. 1. Experimental setup showing a chain of beads with 
sensors SI, S2 and acquisition facilities. The wave that is 
initiated at the left-hand side of this diagram and passes first 
through sensor SI corresponds to LSTW, the opposite wave 
is named RSTW. 



^ = (5/4)2(v/c)4 cos^ [C/(i?V10)] 



(3) 



where c = (2i?)^/''(K/m)^/2. Although this is a trun- 
cated solution of Eq. ([T]), it is well established that 
one hump [— 7r/2 < ^/(i?-\/TO) < 7r/2] of this periodic 
function represents a solitary wave solution [H. Finally, 
approximating the spatial derivative by finite differences, 
the strain in the chain reads ip — 5/{2R). Thus, the 
force at the interface, F ~ K{2R'ip)^^'^ , and the velocity 
V become 
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Our experimental setup consists of a linear chain of 
steel beads. Beads are 100C6 steel roll bearings with den- 
sity p = 7780 kg/m^, Young's modulus Y = 203±4 GPa 
and Poisson ratio v = 0.3. The chain is made of either 

= 26 or = 25 equal beads (radius R=Vi mm) with 
free boundaries, as shown in Fig. [TJ The beads, barely 
touching one another, are aligned on a horizontal Plex- 
iglas track which eliminates transversal degrees of free- 
dom. Two solitary waves are initiated from the impact of 
two small strikers [Rg = 4 mm) located at each ends of 
the chain. Strikers are released from a determined height 
by using an electronic circuit that allows control over the 
time of fall of each striker and consequently, over the syn- 
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chronization of the entrance of each pulse in the chain. 
The pulses are monitored by measuring the dynamic force 
with two piezoelectric transducers {PCB 200B02, sensi- 
tivity 11.24 mV/N and stiffness 1.9 kN//j,m), each being 
inserted inside a bead cut in two parts, as presented in 
a previous work Q. The total mass of an active bead 
and its overall rigidity (led by the contact stiffness with 
neighbors) both matches the characteristics of an original 
bead. The embedded sensors thus allows non-intrusive 
measurements of the force inside the chain. 

In order to register solitary waves interactions and the 
further resulting mechanisms, such as secondary solitary 
waves generation, both sensors are placed in symmetrical 
positions with respect to the middle of the chain (as 
shown in Fig. [1]). These positions are chosen within two 
compromises namely, sensors must be far enough from 
the free boundaries, in order to avoid registering the 
transient formation of the incident solitary waves, and 
they must be separated by a distance larger than the 
typical width of a pulse (about five beads), in order to 
distinguish the secondary solitary waves well separated 
from the main pulses. 



EXPERIMENTAL OBSERVATIONS 

Experimentally, care is taken to establish the exact 
zone where pulses interact. Sources of errors are mainly 
due to amplitude differences of the left-hand side trav- 
eling wave (LSTW) and right-hand side traveling wave 
(RSTW) and the presence of small distance gaps be- 
tween beads which may introduce path differences be- 
tween waves and lead to an undesirable shift in the inter- 
action zone. The later is controlled by softly and carefully 
compressing the chain by hand in order to eliminate any 
distance gap between beads before any acquisition. For 
each releasing of the striker beads, which are dropped 
from the same height, a large set of data that consists in 
the force measurements of incident and passed-through 
resulting waves, is recorded from both sensors. When 
the incident waves registered by the sensors had ampli- 
tude difference no greater than 5%, it is considered that 
solitary waves have crossed at the middle of the chain 
(or, at most, at a distance corresponding to one bead's 
radii from the middle). Data for which this condition is 
not fulfilled are not retained for further analysis. When 
this condition is fulfilled, the reference solitary wave is 
obtained by initiating a single solitary wave with same 
incident amplitude, by dropping a single striker from the 
same height. The reference wave travels along the chain 
from left to right without any interaction. The reference 
data set contains the passing of the reference solitary 
wave through both sensors, from which we obtain the 
unperturbed wave's time of flight, Tg. 
In Fig. [21 both waves are identified as they travel along 



the chain in opposite directions. Indeed, after the colli- 
sion, both waves (RSTW and LSTW) are registered in 
opposite sensors and eventually leave the chain by pro- 
ducing a small gap at the end beads. A time shift appears 
between the resulting waves and the reference wave, and 
we observe that the later arrives to the sensor with a time 
delay A, although the reference wave has a larger ampli- 
tude (thus greater speed) than the resulting waves after 
the collision. The measurement of A is independent of 
the position of the sensors given that it is obtained from 
the subtraction of the flight times of the perturbed and 
unperturbed waves, both being measured at the same 
location. Still, we observed the formation of secondary 
solitary waves (shown in plots (C) and (D) of Fig. [2|) in 
chains with either odd or even number of beads. These 
secondary solitary waves are captured with smaller am- 
plitude in the numerical simulation as shown in Fig. [S] 
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FIG. 2. Experimental result showing the time delay for odd 
and even chains, left column corresponds to a chain of 25 
beads and right column to 26 beads. Plots (A) and (B) show 
the force felt by sensors in positions n — 9 and n = 17 for odd 
{N — 25) chain, while plots (C) and (D) correspond to force 
in n = 10 and n = 17 for the even chain {N = 26). In all 
cases, the solitary wave traveling from left to right (LSTW) 
in our setup is plotted in solid line, dashed line indicates that 
the wave travels from right to left (RSTW), a reference wave 
traveling from left to right is shown with dots. (Notice that 
maximum forces of SW are different for odd and even chains) 



DISCUSSION 

As shown in Fig. HI the time delay between a collided 
wave and the reference is proportional to Fm^^^. This 
observation indicates that the effective interaction 
between solitary waves is a function proportional to 
the duration of solitary waves, r. By considering a 
reference frame over one of the incident solitary waves, 
the interaction itself can be seen as if this wave passes 
over a localized compression zone inside the chain. This 
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FIG. 3. Simulation results (sensors not included). Same nota- 
tion used for Fig. [2] It is observed that viscoelastic dissipation 
does not completely accounts for the decay in maximum force. 
Also, SSW appear in plots (C) and (D) with amplitudes much 
smaller than observed in experiments for the same incident 
waves. 



local change in static load produces an increase in the 
amplitude of the wave and thus an increase of its velocity. 
However, the numeric realization of this idea requires 
the introduction of an estimation of a characteristic 
interaction distance, which might lie in the range of 
the typical width of a pulse, 2R^/T0. Another more 
consistent alternative view is to consider the incident 
waves as two beads approaching with velocity 2Vm one 
to another Q. As before, Vm is the maximum velocity 
of the beads during the solitary wave propagation. 
Since the collision time is a function of relative velocity 
between beads, the time delay is simply given by the 
collision time of a single bead impinging a static one with 
speed Vrm minus the collision time for a frontal colli- 
sion where beads have relative velocity 21/^ 23|. This is. 



numerical calculations. However, in Fig. |4] we observe 
a small difference between experiments and simulations 
for the measured time shift. This difference corresponds 
to the increase in wave duration due to friction dissipa- 
tion between the plexiglas track and beads. This effect 
implies that for a given incident force, pulse duration is 
slightly longer in experiments than in simulations, the 
actual interaction between waves in experiments is thus 
longer than in simulations. 

Fig. [5] shows the SSW amplitudes as function of in- 
cident wave amplitude. The amplitude of SSW is pro- 
portional to the amplitude of incident waves. However, 
secondary solitary waves are stronger when the collision 
occurs at the middle contact in chains with even number 
of beads. This behavior has been predicted by Sen and 
coworkers [l^ who showed that SSW are generated by a 
subsequent release of a little unbalance of potential en- 
ergy at the region where interaction takes place. For the 
case of solitary waves interacting at the contact between 
two beads, the energy unbalance is larger. 

Although our numerical simulations correctly predict 
the existence of these waves, in qualitative agreement 
with results in Ref. [l^, experiments show that SSW 
amplitudes are significantly larger than predicted. This 
difference might have several sources, one of these is the 
energy dissipation that introduces a new time scale in 
the system. However, a first mechanism of dissipation is 
due to viscoelastic behavior of the bead material which 
is included in simulation and does not explain the large 
amplitude of SSW. The second dissipation mechanism is 
due to friction between the plexiglas track and the bead, 
the contribution of friction also affects the duration 
and amplitude of the wave. However, as it was shown 
in previous work this effect is of the same order 
than viscous dissipation in our setup. Thus, neither 
viscous nor frictional dissipation can account for the 
large amplitude SSW observed in experiments. 



2. 



.3\ 2/5 



94 /57rpi?3\ 



Vr 

0.7612 



1 - 



21/5 



V 6 

Cbal X F~^/^ 



en 



4x1/3 



p-1/6 



(6) 



Here the pre-factor, Cbai = 49.74[^sN^/^], is obtained 
explicitly from Eq. ([5]) for the ballistic collision approxi- 
mation. This results is constructed under the assumption 
of elastic interaction including the viscoelastic dissipation 
effect. It is worth mentioning that viscoelastic dissipation 
does not contribute to the time shift because such dissi- 
pation increases the interaction duration of two colliding 
beads in a factor rj independently of the intensity of the 
impacting velocity. Thus, this effect cancels out when es- 
timating the time shift. This result is well verified in our 



It is worth mentioning that recent numerical achieve- 
ments revealed that the presence of a viscous contribu- 



tion, proportional to the relative beads velocity [251, ]2. 



(as it would be for instance in presence of a purely viscous 



fluid between beads [15|, [16| ) or in the form of heuristic 
term proportional to some power of the absolute value of 
the relative beads velocity [27[ , in addition to the nonlin- 
ear Hertz potential, may induces spontaneous secondary 
waves as a result of a slow relaxation of the viscous stress. 
However, the typical time scale of the pulses 25, H] no- 
ticeably differs from the leading pulse, which were not 
observed in our experiments. The model presented in 27 1 
does not conserve momentum by construction and might 
reveal in turn the existence of a complex mechanism in- 
volving losses of momentum through a frictional coupling 
between the beads and external walls or supports. 

In fact, we suspect the difference between the observed 
SSW amplitude and the numerical prediction to be due 
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to small rotation of the beads [28[ induced by the fric- 
tion between the beads and the track. The friction and 
the rotational degrees of freedom are not included in our 
model, despite one can expect that when a single solitary 
wave propagates in the chain of beads, the friction on the 
track induces a small thwarted rotation of each bead. In 
the case where two solitary waves traveling in opposite 
directions interact in one contact, both adjacent beads 
partially roll in over each other, which ultimately leads 
to an excess of potential energy stored into compressional 
deformation. In the other case, when the interaction oc- 
curs at the center of one bead, both adjacent beads roll 
over a static one, which leads to shear deformation of the 
contacts of the static bead. In both cases, the elastic en- 
ergy stored through these additional mechanisms might 
be released afterward and induces SSWs stronger than 
from the frictionless prediction. 
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FIG. 5. Amplitude of the SSW as function of incident wave 
amplitude. Empty markers correspond to experimental data 
while solid markers show simulation results. 
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FIG. 4. Delay between reference and crossing solitary waves, 
A, as a function of maximum incident wave amplitude, Fm- 
Empty bullets and crosses represent experimental data for 
iVo = 25 and A'^e = 26, chains. Solid line represents the sim- 
ulation result and dashed line shows the fit of experimental 
data using A = Ce^p^'™^''^ with Cexp = 45.9 ± 0.5[^sN^/^] 
being in satisfactory agreement with the ballistic approxima- 
tion prediction Chai = 49.74[/isA^''''] (see Eq.[n}. 



CONCLUSIONS 

In conclusion, we have been able to measure the phase 
shift produced by the frontal collision of two equal soli- 
tary waves in one dimensional Hertzian chains. We com- 
pared our experimental observations to the numerical so- 
lution of Eq. (HI). Moreover, by using a well known result 
from linear elasticity theory, we were able to obtain a 
quantitative agreement to account for the time duration 
of the interaction between two colliding spheres with rel- 
ative velocity 2Vm and contrast it to the unperturbed 
propagation where spheres collide with relative velocity 
Vm- Secondary solitary waves have been detected for 
chains with odd and even number of beads. The am- 



plitude difference between experimental SSWs and cal- 
culated ones has been interpreted as consequence of an 
additional potential energy storage due to of solid friction 
at the contact beads. 
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